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Using the formahsm developed in [l], 01 we carry out the first complete calculation of kinematic 
power corrections to the helicity amplitudes of deeply-virtual Compton scattering to the twist-four 
accuracy for a study case of a (pseudo)scalar target. Our main result is that both finite-f, ~ t/Q^, 
and target mass, ~ ni^/Q^, twist-four kinematic power corrections turn out to be factorizable, at 
least to the leading order in the strong coupling. The structure of these corrections is discussed and 
a short model study of their numerical impact is presented. 
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I. INTRODUCTION 

It is generally accepted that studies of hard exclusive 
scattering processes may allow one to access a three- 
dimensional picture of the proton in longitudinal and 
transverse plane Q , which is encoded in generalized par- 
ton distributions (GPDs) [1 [1]. One of the principal 
reactions in this context is Compton scattering with one 
real and one highly-virtual photon (DVCS) which has 
received a lot of attention. The corresponding measure- 
ments are planned, e.g., at the future 12 GeV facility at 
the Jefferson Laboratory. 

The QCD description of DVCS is based on the opera- 
tor product expansion (OPE) of the time-ordered prod- 
uct of two electromagnetic currents where GPDs come 
into play as operator matrix elements. In order to access 
the transverse proton structure one is interested in the 
dependence of the amplitude on the momentum transfer 
to the target t = {p' — pY in a sufficiently broad range. 
Since the available photon virtualities are not very 
large, corrections of the type ^ t/Q'^ which are formally 
twist-four effects, can have significant impact on the data 
analysis and should be taken into account. In addition, 
target mass corrections of the type ~ m? /Q^ can be sig- 
nificant at least in certain kinematic regions and have to 
be included as well. 

Such effects are usually referred to as kinematic power 
corrections since they can be expressed in terms of 
leading-twist parton distributions and do not involve 
"genuine" nonperturbative effects due to quark-gluon 
correlations. A well known example is provided by Nacht- 
mann corrections [0| to the structure functions in deep- 
inelastic lepton-nucleon scattering (DIS). On a techni- 
cal level, these corrections arise because of the subtrac- 
tions that are needed to form traceless operators. The 
Nachtmann power corrections have been studied in full 
detail (see e.g. Q) and are routinely taken into account 
in the analysis of DIS data. Another classical result 
is due to Wandzura and Wilczek [8,] who have shown 



that the twist-three structure function g2{x, Q^) for mas- 
sive targets with spin 1 /2 receives a contribution related 
to the leading- twist structure function gi(x,Q^). The 
Wandzura- Wilczek relation is a consequence of Lorentz 
invariance and can be understood as spin rotation in the 
target rest frame [13] ■ 

The importance of taking into account kinematic 
power corrections to DVCS was acknowledged by many 
authors 0, [Tl| - [20j . This task is much more complicated, 
however, because in addition to Nachtmann-type contri- 
butions related to subtraction of traces in the leading- 
twist operators one must take into account contributions 
of twist-four operators that are related to the leading- 
twist ones by total derivatives. We call these operators 
descendants of the leading twist ones, because their ma- 
trix elements do not involve any new nonperturbative 
parameters. Schematically 



O2 - d^'O,. 



(1) 



where Ou 



are the usual leading-twist operators. The 



problem arises because matrix elements of the operator 
O2 on free quarks vanish [2l[. Thus in order to find 
its leading-order coefficient function in the OPE of e.g. 
two electromagnetic currents one is forced to consider 
either more complicated (quark-antiquark-gluon) matrix 
elements, or stay with the quark- antiquark ones but go 
over to the next-to-leading order. In both cases the real 
difficulty is not the calculation of the relevant Feynman 
diagrams, but the necessity to separate the contribution 
of interest from the contributions of "genuine" twist- 
four operators the number of which increases with the 
spin. This problem was solved in Refs. Jj, t2i using her- 
miticity of the evolution equations for the so-called non- 
quasipartonic twist-four operators with respect to a cer- 
tain conformal scalar product (23j . Hermiticity implies 
that the coefficient functions of multiplicatively renor- 
malizable twist-four operators are mutually orthogonal 
with a proper weight function. Since "kinematic" twist- 
four operators are multiplicatively renormalizable by con- 
struction, using this property allows one to separate the 
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kinematic part from an arbitrary twist-four operator in 
QCD 2] and calculate the coefficient functions for all de- 
scendants of the leading twist operators to the product 
of two electromagnetic currents to the required twist-four 
accuracy [EH- 

The results of Refs.[l|, Q amount to a complete cal- 
culation of kinematic power corrections to twist-four ac- 
curacy in two-photon processes at the level of the op- 
erator product expansion (to the leading order in the 
strong coupling). A question remains, however, whether 
the OPE is applicable to the study of exclusive reactions 
with one real photon, e.g. DVCS. In other language, 
the question is whether kinematic corrections ^ 
to DVCS can be taken into account consistently in the 
framework of collinear factorization. The answer is not 
obvious and, in fact, there are many arguments suggest- 
ing that collinear factorization does not hold in DVCS 
beyond the leading twist. The hope is that the collinear 
factorization framework may nevertheless be valid for a 
subset of "kinematic" power corrections defined as the 
contributions of the descendants of the leading-twist op- 
erators to the OPE, because such terms are intertwined 
with the factorizable leading-twist contributions by the 
electromagnetic gauge and, more importantly, transla- 
tion invariance. In this work we verify this conjec- 
ture by explicit calculation on the simplest example of 
a (pseudo)scalar target. To this end we derive explicit 
expressions for the helicity amplitudes to twist-four ac- 
curacy including all finite-t and target-mass corrections, 
which turn out to be remarkably simple. An inspection 
shows that these amplitudes do not contain stronger sin- 
gularities than those present already in the leading-twist 
amplitudes and, hence, collinear factorization is not en- 
dangered (at least to the accuracy of our calculation, i.e. 
in the leading order in the strong coupling). This result 
is very encouraging since taking into account ~ t/Q^ 
and ~ vr? IQ"^ power corrections removes one important 
source of uncertainties in the theory predictions for inter- 
mediate momentum transfers, ~ 5 — 15 GeV^, which 
is the most interesting range in view of the planned ex- 
periments. 

The presentation is organized as follows. Sect. 2 is 
introductory. Here we define the kinematic variables, ex- 
plain our notation and quote the necessary portion of 
the results of Refs. [l|, that are employed in the fur- 
ther analysis. In Appendix A we explain why the results 
of [l|, Q that have been derived for fiavor-nonsinglet op- 
erators are valid for flavor-singlet contributions as well, 
without any modification. The calculation of helicity am- 
plitudes for the DVCS on a scalar (or, equivalently, pseu- 
doscalar) target is presented in Sect. 3. In this derivation 
we pay special attention to the restoration of translation 
invariance in the sum of all twists, leaving some details to 
Appendix B. The final expressions are collected in Sect. 4. 
In Sect. 5 we present some model estimates of the mag- 
nitude of kinematic power corrections and summarize. 



II. GENERAL FORMALISM 

For definiteness we will consider DVCS on a pion tar- 
get. The hadronic part of the amplitude is determined 
by the matrix element of the time-ordered product of two 
electromagnetic currents 

j7(^)=9>)7^Q9(^), (2) 

where g = {u, d} is the quark field and Q is the diagonal 
matrix of quark charges Q — e diag{e„, Cd}, e = v47ra. 
In the most general form 

y'd4xy'd4ye-'-+^''^(^^(p')|T{j^'"(a;)jr(2/)}k''(p)) = 
= -z(27r)4j(p + g-p'-g')<.(9,'z',P)- (3) 

Here a, h are isospin indices, are the momenta of the 
initial and final state pion, and g, are the virtual and 
real photon momenta, respectively: 

<Z'2 = , q^ = -Q" . (4) 

Making use of translation invariance one can get rid of 
one integration and define the DVCS amplitude by 
a simpler expression 

At, ^iJd^xe-^^^7r\p')\T{j;'^iz,x)jriz2x)}\7r'^{p)), 

(5) 

where zi , Z2 are real numbers subject to the constraint 
zi — Z2 = 1 and we use a notation 

r = ziq- Z2q' ■ (6) 

Since At does not depend on zi , Z2 one can fix their val- 
ues in some way (e.g. set zi — 1, Z2 = or zi = —Z2 = 
1/2) from the very beginning. It is instructive, however, 
to keep these parameters arbitrary throughout the calcu- 
lation for the following reason. In order to calculate At 
including terms ~ 1 /Q^ we have to employ the OPE for 
the product of currents 

T^.{zuZ2) = r{j--(zia;)jr(^2x)} (7) 
to the twist-four accuracy 

T,u = t(*=') + t(*=^) + + . . . . (8) 

The fact that At does not depend on zi , Z2 (transla- 
tion invariance) is a consequence of the transformation 
property of the T-product: 

T^.izi + a,Z2 + a)^ e»(P^) T^.(zi, ^2) e-^-^^^") , (9) 

where P is the usual momentum operator 
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Equivalently, we can write this transformation law in the 
differential form 

(a.i +9^,)t^.(zi,^2) = [(iPa;),T^,(2i,02)]. (11) 

Here and below [* , *] stands for a commutator. 

The crucial point is that these relations do not hold 
for each term in the twist expansion separately, but 
only in the sum of all twists. This implies that the ex- 
pression for the DVCS amplitude obtained from the con- 
tribution of the twist-two operators alone is not transla- 
tion invariant to the 0(1/ Q'^) accuracy, i.e. it depends 
on the choice of the positions of the currents. Trans- 
lation invariance is restored in the sum with contribu- 
tions of higher-twist operators that are related to the 
leading twist by adding total derivatives. These opera- 
tors were constructed explicitly in Refs. [l^ [l^ [13] for 
twist-three and in Refs. [H, Q for twist-four. This means, 
first, that estimates of kinematic corrections based on 
the Nachtmann-type contributions of leading-twist oper- 
ators alone do not have any physical significance and can 
be misleading. Second, translation invariance (indepen- 
dence on 2:1,^2) of the final answer (in the sum of all 
twists) provides an important check of the calculations. 

A similar observation which is thoroughly discussed in 
the literature (see e.g. [O, [H, Ell ) > concerns the electro- 
magnetic gauge invariance of the amplitude A^u, which 
implies the Ward identities 



or, in the operator form 



(12) 



9'^T^,(zi,Z2) - Z2[iP'^,r^.(zi,Z2)] , 

d''T^u{zi,Z2) = Zl[^P^T,,,(zl,Z2)] , (13) 

where — d/dx^. Similar to the above, the Ward iden- 
tities only hold for the sum of all twists in the OPE but 
are violated for each twist separately. It was checked in 
Ref. [2] that Eqs. (|13l) are indeed satisfied in the sum 
of twist-two, twist-three and twist-four terms to the re- 
quired accuracy on the operator level, which guarantees 
that Eq. ([T^ is automatically valid for scattering ampli- 
tudes from arbitrary targets. This constraint is eventu- 
ally built up in the construction of gauge-invariant, e.g. 
helicity, amplitudes (see below), whereas the translation 
invariance condition applies to each helicity amplitude 
separately. 



A. Kinematics 

We use the two photon momenta, q and q' , to define a 
longitudinal plane spanned by two light-like vectors 



where 



n = q' , h= -q + {l-T)q' , 



(14) 
(15) 



For this choice the momentum transfer to the target 

A=p'-p = q-q', (16) 

is purely longitudinal, which is convenient for calculation, 
and the target (pion) momenta have a nonzero transverse 
component: 

A = —h — rn , 
q — n{l ~ t) ~ n , 

P = l{p + p') = l{h-Tn)+P^. (17) 
The skewedness parameter f is defined as 



P+-P+ 
P+ + v'+ ' 



(18) 



where a+ = (an). The scalar product of the two light- 
like vectors n, n in our normalization is of the order of 



(nn) 



whereas 



IP, 



2(1 -r)' 
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(19) 



0(Q°). (20) 



Note that t < and the condition \P\^ > translates 
to the lower bound \t\ > |t„ii„| = 4m2^V(l - ^^). cf. 0. 



B. Helicity amplitudes 

Construction of helicity amplitudes is simplified con- 
siderably in the spinor formalism. Our conventions and 
notation follow Refs. 0, [12, [2^. In this approach each 
covariant four- vector a:^ is mapped to a hermitian 2x2 
matrix Xrvn- 



(21) 



where cr'' = (1, <?), 5-^ = (1, — ct) and a are the usual 
Pauli matrices. We accept the following rule for raising 
and lowering of spinor indices 



u 
u° 



1 I . C -11^ 



(22) 



where the antisymmetric Levi-Civita tensor e is defined 
as follows 



ei2 



,12 



=12 



1. 



(23) 



For this definition 
S'j and (e"'')* = e^'^ . 
spinors is defined as: 

(uv) — U°'Va — —UaV" , [uv) — UaV" 



SP and e"^ = -e/ 



An invariant product of Weyl 
-u"va. (24) 
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More details and some useful identities can be found 
in Hi-^]. 

An arbitrary like-like four-vector in the spinor repre- 
sentation, Oq-q, = 0, can be parameterized by a Weyl 
spinor such that aaa = £,a£,6n where = SX- particu- 
lar we introduce two auxiliary spinors, A and /i associated 
with the light-like vectors n and n defined above: 



(25) 



Thus the two-photon amplitude Afn, is parameterized 
in the most general case by nine scalar functions A'^^''^^] 
three of them correspond to a longitudinal photon in the 
final state and, therefore, do not contribute to the cross- 
section. 

In order to bring this expression to a more familiar 
form we define three photon polarization vectors, ~ 
= 5<"£od ) as follows 



For example one can choose 

The scalar product (nn) can be written in the form 

1{nn) EE 2(n^n^) ^ n^an'^" = i^^>^){>^P■)■ (27) 

One of the advantages of this formalism is that two aux- 
iliary spinors are sufficient to define a basis in the whole 
four-dimensional space. Indeed, the basis vectors in the 
transverse plane (orthogonal to n, n) can be chosen as 
/IqAq and Aq/2q. An arbitrary four- vector can be ex- 
panded in this basis as 

2{nh)xaa — (28) 

= X^^fiafia + X XaXa + X — ^-Xa^a + — MaAd, 

where 

X-+ — fl°'XadX°' , X = fl°'XaaP'°' (29) 

so that and x correspond to the "plus" and 

"minus" coordinates in the usual sense, respectively, 

whereas x^ and x |- are the two (holomorphic and 

anti-holomorphic) complex coordinates in the transverse 
plane. 

Using the spinor formalism it becomes straightforward 
to write down a general parametrization for the DVCS 
amplitude 



(30) 



that takes into account the Ward identities p^ . The 
constraint q'^A^,, = becomes A^A^^^^^^ = 0, and its 
general solution is 



(2) 



A 



(3) 



(31) 



The first two terms in this expression correspond to the 
contributions of a transversely polarized real photon in 
the final state, and the last term to a (unphysical) lon- 
gitudinally polarized photon. Making use of the second 
relation in Eq. (|17l) one can easily resolve the other con- 
straint, q^A^i, = 0, that results in 



Ai'i = A„/2AA(^'i)+/i„Ad^('''^ 

+ iHaP-a + (1 - T)AaAd)^(' 



3) 



(32) 



_0 _ Hafia + (1 



\/{nn) ' 



^2(1 -T)(nn) 



(33) 



The polarization vectors are normalized as (e^e '') = 
— 1, (e^)^ = 1 and are orthogonal to the photon mo- 



menta, q^e^'^ 



0. Note that (e+)* 



Thus we can write, finally. 



A — F+F^A'^ 



,0^+ /io-_L=-+=-+_4+- 



' A 



-£ £ 



^"e-A"+ 



tA<;!\ (34) 



The amplitudes A^^ (A^^) describe virtual Compton 
scattering of transversely polarized photons of positive 
or negative helicity without (with) a helicity flip. The 
amplitudes A'^'^ correspond to the longitudinally polar- 
ized virtual photon in the initial state. The amplitude 
A'^'^^ is of no interest as it does not contribute to any 
physical observable. 

Further constraints follow from parity conservation. 
One obtains 



A++ 
A^^ 



A— = A'^"^ 
{e^P^fA^^' 



(35) 



To explain the last two equations, note that the light-cone 
vectors vectors n, n, Eq. ([25]), are invariant under the 
redefinitions of the auxiliary spinors by arbitrary phase 
factors 



A 



s'-^iA and 



(36) 



The amplitude A^u must stay invariant under such rota- 
tions as well. Taking into account that with our choice 
of kinematics the only vector at hand which has nonzero 
transverse components (and, therefore, a nonzero scalar 
product with the polarization vectors e^) is P'^, one eas- 
ily arrives at the representation in Eq. ([35]). 

We end up with three independent helicity amplitudes, 
A^^\ A^^\ A^'^\ (or, equivalently, A++ , A°+ , A-+) and 
will calculate them in what follows. 

The amplitude A^v in vector notation can be written 
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in terms of the hclicity amplitudes in the following way: 



where daa = cr'^ . 9„ and 



1 



(37) 



where 



im') 



quq^ 



iqq') 



l\2 



a '13 



(38) 



The last term ~ ql in ((37| does not contribute to physical 
observables. 



C. Operator product expansion 

The OPE for the product of two electromagnetic cur- 
rents, Eq. ([7]), to the twist- four accuracy taking into ac- 
count the leading-twist operator and its descendants ob- 
tained by adding total derivatives is available from [l|, [3] . 
In this Section we collect the necessary expressions. 

To be precise, the twist-four contributions to the OPE 
have been obtained in P, @) for flavor-nonsinglet opera- 
tors only. We are able to show, however, that the same 
expressions are valid (at tree level) for flavor-singlet op- 
erators as well. This generalization is explained in Ap- 
pendix A. 

The tensor 2^jj^(j^(-2i, 22) can be written to the twist- 
four accuracy in the following form: 



a/3al3 



2 r 



2 

Xl3a^ap{z2, Zi) + X^^Xpa AA(zi, Z2) 

C(zi, 22) - Xaadpp£{z2, Zl)j | , 

(39) 



where zi2 — zi — Z2. (We will set zi2 — >■ 1 in the final 
expressions) . 

The operators AA and C are twist-four, whereas the 
expansion of *8 starts from the leading twist: 



*R ■ — -t- -I- 



(40) 



The twist-two result is well known and can be cast in the 
form 



1 

*8ad'(2l,22) = J duO\ 



l(uZi,UZ2), (41) 



O^f (zi,Z2) = [nO++](zi,Z2) . (42) 
Here 11 is the (leading- twist) projector defined as 



Aq 



A=0 



(43) 



(see Ref. [1] for details) and 0++{zi, Z2) is a certain com- 
bination of vector and axial- vector light-ray operators 



0++(zi,Z2) = ■^[Ov{zi,Z2) ~ Oviz2,Zi) 



Oa{zi,Z2)-Oa{z2,zi)\ (44) 



which are defined as 



Ov{zi,Z2) = q{zin)Cf- ijiq(z2n) , 
Oa{zi,Z2) = g(zin)Q^ ^75g(z2?^) 



(45) 



We remind that Q is the matrix of quark electric charges, 
Q2 = e2[5/18 + 1/6t3], cf. ©. In the pion matrix ele- 
ment of the operator 0++ only the vector isoscalar con- 
tribution survives: The pion matrix element of the axial- 
vector operator Oa vanishes identically, whereas the ma- 
trix element of the isovector operator 0^^{zi , Z2) is sym- 
metric under permutation of the arguments and drops 
out from the expression in Eq. (|44p . Thus for the case at 
hand 



0++(zi,Z2) = >c u{zin)'^u{z2n) + d{zin)'^d{z2n) 



(46) 



where 



5e^ 
18 



(47) 



Because of this structure, the expression for scattering 
amplitudes from a pion are identically the same as for 
a scalar target apart from a trivial Kronecker symbol in 
isospin indices, A'f^'l = S'^^A,, 



The twist-three contribution to the OPE is contained 
in S*„=3(zi,Z2): 
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udu 



Z2 ^12 



(48) 



Note that the last term in this expression is by itself 
twist-four; its role is to subtract twist-four contributions 
from the first two terms so that the result is a purely 
twist-three operator. 

The twist-four contribution to *Bqq(zi, Z2) can be writ- 
ten as 



*BaA(-^l:^2) = x'da,&M{zi,Z2), 



(49) 



SO that the complete twist-four contribution to the OPE 
is determined by three functions: B(zi, Z2), C(zi, Z2) and 
AA(zi,Z2) = A(zi,Z2) — A(z2, Zi). Explicit expressions 
for these functions are given in Ref. [2] in several equiva- 
lent representations. We found the following form to be 
the most convenient for our present purposes: 



(zi,Z2) — — f du < \nU ZiZ2 0i{ziU, Z2U) 

4 Jo I 

+ (z2dz2 - — - \nu Z2d'^ Zi2]TZ{uZi,UZ2) 

V 2:12 ^ 

— (zidzi — \nu zid'^^Z2i)'JZ{uzi,uz2) 

^ Z21 ' 

(50) 



1 f^du 

{zi,Z2) = - / — \ u^{l-u'^+u^lnu)ziZ20i{ziU, Z2U) 



TZ{uZl,UZ2) 



TZ{uzi,uz2) >, (51) 



(1 - U') Z2d, 



^^2 

212 

,2 I „,2 , „,N „_p,2 



(l-u^-f-u^ \x^u)z2dl zx 



Z2 



(1 - u^) 

V ^21 

-I- (1 — Inw) z\d\^Z2\ 



and 



C(zi, Z2) 



- / '^\TI(uz\,uz2) +7?.(uz2,'uzi)] , (52) 
8Jo " 



where 

Tl{zi,Z2) 



Z\2 



dw 



dw' w' — Z2 



Z2 Zl2 Zi2 Zi 



-S+Oi{wM)-{So-l)02{w,w') 



'R{zi,Z2) = Z12 



dw 

Z\2 



dw' zi ~ w 

Z12 W ~ Z2 



-S+Oi{w,w')-{So-l)02{w,w') 



.(53) 



Here S+, So are differential operators acting on the light- 
cone coordinates 



S+ = w'^dr, + 2w + w''^dyj' + 2w' , 
So — wdw + w'dw' + 2 



(54) 



and 01^2 are nonlocal (light-ray) twist-four operators de- 
fined in terms of the leading- twist operator 0^(21, Z2) 
as follows: 



Oi{w,w') 

02{W,W') 



0'-^{w,w') 



(55) 



The matrix elements of these operators are expressed in 
terms of the leading-twist matrix elements (i.e. GPDs) 
as 

{p'\0,{w,w')\p) = -A2(p'|0^tK^')|p), 

d 



{p'\02{w,w')\p) 



As mentioned above, the pion matrix element of 
the axial-vector operator vanishes, so that effectively 
0++ (2^1, -22) is antisymmetric under permutation of the 
arguments for the problem under consideration. As the 
result the operators 01,2(^1,22) are antisymmetric in 
zi, Z2 as well, whereas 7?.(zi, Z2) = 'R-{z2, 21 )■ 

The reason why this representation for the twist- 
four contributions turns out to be more convenient 
as compared to more explicit expressions in terms of 
^1,2 (see [2]) is that the 7?,-operators are themselves 
translation-invariant to the twist-four accuracy: 



{dzi +9^2)7^(^l,Z2) = ^[(Px),7^(zl,Z2)], 



(57) 



see Appendix B. This property simplifies the structure of 
rather intricate cancellations that are otherwise necessary 
to restore translation invariance of the final results. 
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D. Pion GPDs 

The matrix element of the operator 0+_|_(zi, Z2) sand- 
wiched between pions with different momenta can be de- 
fined in terms of the (isoscalar) pion GPD 25-27] 

{n\p')\0++izin,Z2n)\n^p)) = 



In this form, which turns out to be the most convenient 
for our analysis, the pion matrix element is parameterized 
in terms of a single DD 

a, t) = dpfiP, a, t) + a, t) , (64) 

which is related to the pion GPD H{x, ^, t) as 



2P+6-''k da; e-»^+[-i(«--)+-^(-+«)li/(a;,e, 0,(58) d,H{x,C,t) = f dad/S 6{x - (3 - ^a) m3,a,t) . (65) 



where m: is defined in Eq. (1471) . 

An alternative parametrization of the same matrix el- 
ement is via the double distributions (DDs) [2^, [28| . 
For the isoscalar case one usually introduces two pion 
DDs [29II, /(/3, a, t) and g(/3, a, t), so that 

{^\p')\0++{z,n,Z2n)\^''{p)) = 

1 1-1/31 

= le 

-1 1/31-1 



''^"x(5''^ I d(3 I da 6-^(^^1-2") 



X 



2(Pn)/(/3,a,t)-(An)5(/3,a,t) , (59) 



where 

4, 



-ziA + iz2-zi) /3P--(a + l)A 



(60) 



The isoscalar pion GPD is given in terms of the DDs by 
the following expression: 



H{x, = / dadf3 S{x-P-^a) f{P,a, t)H9{P, «, t) 



(61) 

which can be obtained by inserting 1 — J dx S{x — /3 — S^a) 
under the integral and comparing the definitions. 

Charge conjugation and time invariance impose the fol- 
lowing symmetry properties of the DDs [soj : 

f{/3, a, t) = /(/3, -a, t) , .g(/3, a, t) - -.g(/3, -a, t) , 
/(/3, a, i) = -/(-/3, a, t) , g(/3, a, t) = .g(-/3, a, t) . 

(62) 

The description in terms of two DDs is, however, 
redundant (see a detailed discussion in [s^l), because 
the expression for the matrix element is invariant under 
their redefinitions. This freedom can be used to reduce 
g{f3,a,t) to some minimal form (the "D-term" [1^), but 
it is more convenient (s 
definition in Eq. ((59|) as 



it is more convenient (see e.g. [30|, |31|) to rewrite the 



(the 

m 



{Tr'{p')\0++{z,n,Z2n)\7T-{p)) 
= —xS"'' [ dpda 



Z21 
2i 



fdp + gda 



Z12 



dpf + dag 



, (63) 



where for the last representation we assumed that /, g 
vanish at the boundaries. 



As a consequence of the symmetry properties (15^ the 
function $(/3,a,i) is even under reflection (/3,a) 
(-/3, -a), i.e. 



$(/3,a,i) = <i>(-/3,-«,i) 



(66) 



We note here that the assumption that the DDs / and 
g vanish on the boundary is made only for convenience 
and can be relaxed without any effect on the final result. 



III. CALCULATION OF HELICITY 
AMPLITUDES 

Before going into details of the calculations we want 
to make some general remarks. Note that definitions of 
the pion GPDs involve matrix elements of nonlocal oper- 
ators at a strictly light-like separation, whereas the OPE 
is written in terms of operators with arbitrary (non-light- 
like) separations which involve the leading-twist projec- 
tors: 

{TT\p')\0'^^{z^x,Z2x)\n%p)) = 

= I{{x,n){TT\p')\0++{zin,Z2n)\TT''{p)) . (67) 

Since the dependence of the matrix element on the quark 
coordinates is entirely through the exponential factor 
g-i(^ziz2")j cf. the leading- twist projector is effec- 

tively applied to this exponent. It is sufficient for our 
purposes to know the two first terms in the light-cone 
expansion [33i] 



n[e 



]{x)=^ 



dvve-'"^^''^ +Oix^). 

(68) 



Thus one obtains 



{p'\OX+{zix,Z2x)\p) = 



^4 



2ix 
Z12 



dfida $(/3, a, t) 



(69) 



Here and below we suppress the isospin indices. Note 
that 



Z\Z2 + 2120(212^" - Zx) 



(70) 
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where we introduced a notation 



F : 



K?— ) 



(71) 



The second term 



in Eq. involves correc- 



tions that are proportional to the momentum transfer 
t = A'^ and target mass squared m^, cf. (PO)) . The ex- 
tra factor is converted to a suppression after 
the Fourier transformation ([5]), so that these contribu- 
tions give rise to power suppressed, rv? jQ^ and i/Q^, 
corrections to the hehcity amplitudes. They correspond 
to the generalization of the Nachtmann corrections for 
DIS to off-forward kinematics and have been discussed 
e.g. in iHIil. 

To the same accuracy one must take into account con- 
tributions of "kinematic" twist-four operators ((55|) . The 
leading-twist projection operator has no effect (to our 
accuracy) in the case of Oi, since the matrix element is 
itself proportional to , but it must be included for O2 
because of the 9/9a;p-derivative, see One obtains 



(p'|Ol(zi,Z2)|p) 
(p'|02(^l,^2)b) 



A^y" dad/3 $e-*(^^i^2^^\ 
dad/3 $[(A4i2je-'(^^i'-2^=") 



2i>f 
Z12 



where 



Z2 - azi2\ 



(72) 



(73) 



The second term, ~ {iAx)i1_^^^, in the expression for 
{p'\02\p) is due to the twist projection. Note that this 
contribution contains terms ^ (through the depen- 
dence on I Pip, cf. Eq. ([70]) ) and, therefore, the calcula- 
tion of both target-mass and finite-t corrections requires 
taking into account twist-four operators. We emphasize 
this point because, naively, "kinematic" twist-four op- 
erators in Eq. ([l} are only relevant for the momentum 
transfer dependence. 

Factorization is proven for the DVCS amplitudes to the 
leading power accuracy only and there are reasons 
to expect that the OPE does not give a complete answer 
beyond the leading twist. The main question that we ad- 
dress in this study is whether a subset of power correc- 
tions vr? IQ^ and t/Q^, which we call kinematic, defined 
by the sum of contributions of the leading-twist opera- 
tors and their descendants in the OPE, is well-defined 
for the particular kinematics of DVCS, where one of the 
photons in on-shell. Our calculation will show that the 
tree-level kinematic corrections to all helicity am- 

plitudes are finite, and their form is consistent with fac- 
torization. Whereas this exercise does not constitute a 
proof, it gives a strong argument in favor of the factor- 
ization conjecture and allows one to take into account 



kinematic corrections, which can be numerically large, 
within the conventional framework. 

In order to have the calculation under control at the 
intermediate steps it is convenient to consider a more 
general process with two virtual photons, in which case 
OPE can be expected to hold to all twists, and take the 
limit (q'Y — at the end. Since the photon momentum 
q' enters the calculation via the combination of momenta 
r = ziq — Z2q' ~ q' + ziA only, see Eqs. ([S]),®, we can 
redefine 



ziA 



r = q 



CA: 



(74) 



and take the limit ^ = zi at the end of calculation. This 
change is equivalent to a redefinition q' — >■ q' + — zi)A 
which implies q'^^Q and effectively regularizes all inte- 
grals. 



A. Helicity flip amplitudes 

We begin with the calculation of single- and double- 
helicity-flip amplitudes 



A0+ _ + Afiu 

In chiral notations 



A- 



4(rtn 



-A- 



A A^ 



(75) 



(76) 



and 



A 



0+ 



1 



1 



4(nn) ^2(1 -t) 
1 1 



2(nn) ^2(1 -r) 



[A. 

A- 



-+ + (1 - t)A+-++\ 
+ , (77) 



(1 - t)A+- 



due 



where we have used that A 1- 

to the Ward identity p2)) . 

In order to find these amplitudes one has to insert the 
OPE ([39]) into Eq. ([5]) and evaluate the corresponding 
integrals. Note that because of their symmetry proper- 
ties under rotations in the transverse plane the helicity- 
flip amplitudes must be proportional to powers of the 
transverse momentum A^^ oc Pi, A ^ cx P^, multi- 
plied by scalar functions that depend on A and P_l at 
least quadratically, cf. Eq. (|35l) . If the calculation is 
done to the accuracy, this dependence is irrelevant 

and can be neglected. As the result, only twist-two and 
twist-three operators, *8*=^ (gl]) and S*^^ (gSl), have to 
be taken into account 



A- 



= A'l 



A 



t=3 



(78) 



and similar for A [-+. Twist- four operators do not con- 
tribute (to this accuracy), which is a major simplification 
compared to the case of the helicity-conserving amplitude 
which will be considered in the next Section. 
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By the same reason it is sufficient to keep the first (triv- 
ial) term in the expansion of the leading-twist projector 
n only, see Eq. (|68l) . In addition, one can safely neglect 

all terms which contain a "minus" derivative: d — 

(pdfl) = 2{nd). Indeed, when applied to the matrix 
element, it produces terms ~ {n,i) ^ (niP), (nA) ^ 
0{t,w?). Taking the matrix element of Eq. (|t5|. 

the action of the momentum operator P'' is effectively re- 
placed by the multiplication by n^: P — >■ A = —h — Tn ~ 
—n. Thus the second line in Eq. p8)) gives rise to contri- 
butions 0{t, m? ) only and can be dropped. It is also easy 
to check that only the terms proportional to Z2 in the first 

line contribute to A i-; all other terms contain d 

and can be omitted. 

One finds after a short calculation 



A 



t=2 



A 



t=3 



4m: I dad/S ^{l3,a) I duu 







X Z2 



dw 



zi2 {(r - 



(79) 



where we used that {^£ziZ2^) — ^Pz^P^x since A does 
not have a transverse component. 

Further, to our accuracy r^^ = 2{nn) and 

(r -I- u^tu^j)^ = —2{nh) [zi — uw ~ uiz^ ~ w)F\ , (80) 

where F is defined in Eq. ((7T|) . Taking the integrals over 
u and w one obtains 



A 



t=2 

h 



/ dad/3 <I>(/3,a)/3 



ln(i^— ie) — In Z\ 



A 



t=3 



dad/3 $(^,a)/3 
'ln(F-ie) 1 



F 

Z2 



Z\ 



F^zi 



F - 1 



Z2 



In zi 



(81) 



Note that the contributions of twist-two and twist-three 
operators are separately not translation invariant, i.e. 
their contribution depends on zi,Z2- This dependence 
cancels, however, in the sum, as expected. We end up 
with a very simple expression (cf. |13l - [l5| ) 



V2P, 



dadp<P{l3,a)l3 ■ 



In [F - ie) 
F-1 



(82) 



which is the final result. 

Calculation of the double-helicity-fiip amplitude A ^ 
is similar. Inserting the OPE (15^1) in Eq. ^ and evalu- 
ating the Fourier integral one finds 



A': 



8m dadf3^{fi,a) I duu 



(r + u4iz2)^ + *e 



A 



t=3 



dadl3^{(3,a){pP^x) / duu-" 





zi{w - Zi) 



Z2{Z2 - W) 



{{r + u£z,^)^+ie)^ iir + ue^z2)^+ie)^ 



dw 



.(83) 



After performing the u and w integrations the results can 
be written as 



A 



t=2 





(nfi) 



dadp $(/3, a) 0^ 



F zi 

\n{F — ie) — — In zi 



A 



t=3 



-++ 



p 

[nh) 



dadl3<t>{l3,a) ^^df 



In (F - ie) 



l-F 



2zi 

F-zi 



\n{F — ie) — In zi 



, (84) 



where dp = d/dF. In this case, again, we observe that 
translation invariance is restored in the sum of twist-two 
and twist-three contributions: 



A-+ = - 



2(?in)2 
X dp 



dad^^{l3,a) ^2 
2F — 1 

In (F - ie] 



F- 1 



(85) 



The both amplitudes, ([5^ and (1551) . can be rewritten 
in terms of the "standard" GPD H{x,^,t). The corre- 
sponding expressions will be given in Sect. IIVI 



B. Helicity conserving amplitude 

Calculation of the helicity-conserving amplitude 

1 



A^+ = F^F+A'^" = 

^ " 4(nn) 



(86) 



to the 0{1/Q^) accuracy presents our main goal. This 
calculation is much more involved because of the contri- 
butions of twist-four operators. We can write 



A++ = A^ 



(87) 



where AfJ^^ presents the main challenge. 

Instead of calculating the A-\ spinor projection of 

the OPE directly, it proves to be convenient to proceed 
as follows. Consider a trace of the DVCS amplitude in 
the Lorentz indices 



1 



4(nn 



(88) 



Taking into account that = A^^ and A | |- = 

thanks to (fT^ . we can express the contribution of interest 
as 



A 



++ 



1 



8(7in) 



A+-+- - A-+- 



(89) 



This representation turns out to be very convenient be- 
cause of strong cancellations of different contributions in 
the trace and in the difference of the two terms with dif- 
ferent polarizations in the square brackets. 
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Taking the trace of the OPE m Eq. p9l) one obtains 
1 



It*' 
2 ^ 



2^,4^3 



12 



2;^A<B^(zi, Z2) + x^AA{zi, Z2) 



x^(x9)AC(zi,Z2) 



(90) 



where 



A*B^(zi,z2) = Q3^(zi,z2)-«B^(z2,^i), 

AC(Z1,Z2) - C(Z1,Z2)-C(Z2,Z1). (91) 

We remind that A and C operators are of twist- 
four, whereas 03^ contains all twists. However, since 
a;^Q3^"^(zi, Z2) — (by construction), only 03*=^ 
*8*^'* operators contribute to the trace which, therefore, 
receives no twist-three contribution at all. 

Making use of the explicit expression in Eq. (pi and 
the identity 

(x9)C'*|!^(mZi,UZ2) = {udu + l)C'++(wZi,WZ2). 

we can simplify the leading twist-two contribution to 



X^A^'-^{Z^,Z2) = OXl(.z^.Z2) . 



(92) 



From now on we tacitly assume taking the matrix 
element over a (pseudo)scalar target in which case 
(p'|C'++(zi, Z2)|p) is antisymmetric in zi,Z2. Similarly, 
for the 7?.-operators we use 7?.(zi,Z2) = ^(z2,zi) which 
should be understood as the relation for the correspond- 
ing matrix elements. 

The twist-four contributions to the trace can be simpli- 
fied as well. Taking into account that x^A*B*;^'*(zi, Z2) = 
a;^(a;9)B(zi, Z2), see Eq. one obtains 

(a;a)AB(zi,Z2)-f AA(zi,Z2) = 



du 



ZiZ2U^Ol(uZi, MZ2) — Z2d2Zi2TZ{ziU, Z2U) 



Zldfz2lTZ{z2U, ziu) 



and 



(xa)AC(zi, Z2) = -- 7^(zl, Z2) - 7?.(z2, zi) 



(93) 



(94) 



Next, consider the expression in the square brackets 
in Eq. (|89p . In this case only the antisymmetric in 
(aa) o (/3/3) part contributes to the answer resulting 
in 



c+_ [S*!!3(zi, Z2) + 'B*!!^ (z2, zi)] 
c_+ [«BV=i(zi, Z2) + 'B*ri(z2, zi)] 



+ x^ 



X — d++ [C(zi, Z2) + C(z2, zi)] 
x++9__[C(zi,z2)-t-C(z2,zi)] \.{95) 



The last term (^ x++d ) gives rise to corrections of 

order 0{t^,m^t, m*) after the Fourier transform and can 
be neglected. In the remaining term x"^ one can replace 
X--d+-\. — )■ A{nn){xd) since the difference is again a cor- 
rection ©(t^, m^i, m*), after which it can be combined 
with a similar term C in the trace. 
Collecting everything we obtain 



■^t=2 



d x e 



A++ - - 



TT^ X^ 

1 f d'^x e"*''^ 



{p'\0'-Hz,,Z2)\p), 



A{nn) 



x+^{p'\^'rl{zu Z2) + sit (^2, ^i)b) 

- X_+(p'|'B*p3(zi,Z2) + ®t'(^2,^l)b)' 



d^x e-*''^ , „ 
— —\P 



du 



ZlZ2u'^Ol{uZl,UZ2) 



- Z2dlzi2TZ{ziU, Z2U) + Zidlz2lTZ{z2U, Ziu) 

-27^(z2,Zl)}|p), (96) 

where we suppressed overall l/z^2 factors on the r.h.s. 
which are irrelevant as we have to set Z12 ^ 1 in the 
final expressions. 

We begin with the twist-two contribution. Using the 
representation for the matrix element in Eq. (j69p and 
taking the Fourier integral one obtains 



■At=2 ^ 2 / dad(3m3,a)\\n{{r + £,,,,)^+ie) 



1 f 1 
^5,, 5,^ / udu -7^ - 



(97) 



Note that the second term is already 0{t, m^) because of 
^1-^22^ that we can simplify the denominator as 



(r + u£ziZ2)^ + ie = —2{nn)[uzi + uF — ie 



(98) 



where u = 1 — u, whereas in the first term we need an 
expansion of the logarithm to the 0{t,m^) accuracy 



ln((r + 4i.J'+«e) = HF - ie) 

, 1 f^'\P± 



2{nh) F-ie 



(99) 



The ellipses stand for the terms that are constant or lin- 
ear in a, P, which can be dropped because 



dadp $(/3, a) [a + ba + c/3] = 



(100) 



due to Eq. and the reflection symmetry ^(/S, a) 
<i?(-/3, -a). Using 



udu 



uzi -\- uF — ie 



Fhi{F~ie) zilnzi 



F-zi 



F-zi 



(101) 
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we obtain 



AJ2 = 2 / dadp<i>{l3,a)Un{F-iO) 



Zl 



2{nn) { ^' F 

[ziZ2 + {F- zi)a -~F{F~1)] 
-F\n{F-ie) zilnzi 



X dp 



F-zi 



F-zi 



(102) 



The calculation of the twist-three contribution follows 
the scheme described in Sect. (IIII Al) . We mention only 
that the second line in Eq. does not contribute to 
the result due to the antisymmetry in zi,Z2- Also, in 
this calculation it is sufficient to take into account the 
leading term in the twist projection operator (j68| only. 
A straightforward calculation yields 



A 



1 



t=3 



[nn) 



dadl3<t>{l3,a) 13 

\n{F - ie) 



F ~1 



, (103) 



where \P±\^ is defined in Eq. (|20l) . 

Calculation of the twist-four contribution is somewhat 
longer. It is given by a sum of four terms, 



A++ _ aW , Ay-^) 



(2) 



t=4 



(104) 



which we numerate in the same order as they appear in 
the last equation in (|96|. 

The first contribution, of the operator Oi, is really 
simple. A short calculation yields 



A^'\ 



dadl3^iP,a) 



X d 



(nn) 

F\n{F — ie) zilnzi 



F-zi 



F~zi 



(105) 



The remaining three contributions involve the 7?,- 
operator defined in Eq. (|53|) . As the first step we evaluate 
the integral 

lH(C,zi,z2) = / ^^(p'|7e(zi,z2)b) . (106) 



We take r = q' + (A and will take the limit C — >■ 21 at 
the end of the calculation. The relevant matrix elements 
of the "kinematic" twist-four operators (^'101(2)1^) are 
written in terms of the leading- twist DD in Eq. (l72t . 
Inserting these expressions in Eq. ([5^ and making use of 
the identity 



{So - 1) , / dTOe-*"(^™"''^) 



/2 

ww' ^— M^,„.,,/a^) 



(107) 



we obtain (note that we do not assume Z12 — 1 here) 

f P'^dw' 
{p'\n{zi,Z2)\p) = -2i dadfi^iP.a) dw 

J Jz2 Jz2 ^12 

w' - zo ( 

X U —S+ + {So - l)(A£w) 



Zl — w 



-{Ax)i 



2 

ww' 



r- -(108) 



Going over to the momentum space one gets after some 
algebra 



9^(C,^l,^2) = 



dad/S $(/3, a) 



(nn) 

,2^171 lArt) , a2. 



+ A^{F-l)FdF + A''a 



/(C,zi,z2),(109) 



where 

I{C,Zl,Z2) = 

'■^^ dw 
Z12 



(110) 



dw' 



W - Z2 



Zl — w' C — w — [w' — w)F — ie 



1 /"^^ dw' w' - Z2 ^ C- Zl - {w' - zi)F - ie 



F - 1 



Zi2 Zl ~ W' 



C — w' — ie 



The second and the third twist-four contributions in 
Eq. (|96|) have similar structure, so that they can be 
added together. These terms involve differential opera- 
tors, Z29|zi2 and zidiZ2i, acting on the field coordinates. 
It is easy to show that 



Z2dizi2l{C,Zi,Z2) 
Zidfz2lI{C,Z2,Zi) 



1 Z2 , C,- zi + Z12F - ie 
m - 



F - 1 Z12 



C~Z2 



1 Zl C - ^2 ~ Z12F - ie 
m ■ 



F - 1 Z12 



Rescaling Zi — > uzi we can put Q = zi and Z12 — 1 and 
perform the u-integration. The result reads 



/I (2+3) 
-^t=4 



1 



2{nn) 



dadl3<^{P,a) 



P'\Pi\dF + A'{F - l)FdF + A'a 



1 



F -I 



Z2 



F , Zl F - 1 , 

In F In Zl 



F-zi 



Z2 F - Zl 



Zl 



1 -F 



22 + F 



In Zl 



1 -F 
22 + F 



ln(l - F) 



(111) 



where F ^ F — ie inside Z2[. ■ ■] and F F + ie inside 
zi[...]. 

Finally, for the last twist-four contribution we need 



I{C,z:2,zi) 



ln(l - F~ie) 



F -1 



Li2(l)-Li2(F-fie) 



(112) 
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where Li2(l) — 7r^/6. Changing variables in the first 
term, a — ^ —a, /3 — > — /?, and using that under this 
transformation _F — > 1 — one obtains 



1 



[nn) 



dadp a) 



l3''\Pl\dF + A2(F - 1)F9f + A^a 

^ ^ln(F-ie) + Li2(i^ + 2e)-Li2(l) . 

(113) 



F- 1 



Note that Li2 (a;) is analytic at a; — > but has singularities 
in the expansion around a; — > 1 , which is the reason that 
we keep the +ie prescription in Li2(F + it). 

Collecting all twist-four terms and adding the contri- 
butions of twist-two and twist-three operators we obtain 
the final result: 



X J dadp <^>{I3, a) l^2log{F -ie) 
'l3^\P±fdF-A\F-a) 



1 

(nn) 
1 
2 



1 



F~ 1 

\n{F-ie) + U2{F+ie) - Li2(l)j I. (114) 



This expression is rather remarkable. First, notice that 
the dependence on zi_2 disappeared so that translation 
invariance is restored. Second, the result is finite: there 
are no infrared divergences which might signal break- 
down of factorization to the accuracy because of 
a real photon in the finite state, and would correspond 
to singularities in the limit C — > zi. Third, singularities 
of the 1/Q^ correction at — >■ and F 1 (which cor- 
respond to a; = ±^ in the GPD notation) are not stronger 
than the (logarithmic, in present variables) singularity of 
the leading twist result. This property (see more details 
in the next section) ensures that non-analyticity of the 
GPD at a; = ±^ does not lead to any problems, there are 
no non-integrable singularities in the momentum fraction 
integral. 



IV. FINAL EXPRESSIONS 

In this Section we present our results for the helicity 
amplitudes in terms of the GPD H{x,£,,t) (|65p. Rewrit- 
ing 



F-a= ^ 



1 - F 



we left with integrals of the following structure: 

/i = / dad/3^{P,a)Y{F) , 
h = I dadl3<^{f3,a)l3Y{F), 



h = J dadl3<Pil3,a)(3'dFY{F), (115) 
where Y{F) is a certain function of the F- variable defined 

in Eq. ([7T|l . Inserting 1 — dx 6{x — (3 — a£,) under the 
integral one obtains after a short calculation 



h = 



dxY 



(^)a.iI(x,C,i), (116) 



2e 



h = 



j'^ dxY^^) {xd. + ^di)H{x, t) 

J\xy[^)h{x,U), (117) 



and 



/3 - -^^i_f^y{^ 



id,x+^d^)ixd, + ^d^)Hix,^,t) 

= -2edl l'^dxY[^)H{x,U), (118) 

where we used the identity 

l3S{x-(3- aO = {xd^ + ^d^)e{x - /3 - aO (119) 

to arrive at the last two expressions. 

Using these integrals it becomes straightforward to 
rewrite the results of the previous Section in terms of 
H{x,^,t). We obtain the helicity-flip amplitudes (in the 
notation of Eq. p5p 

H{x,tt)^ fx + ^ 



2e 



Q J x-£, 

^<..-,.a|/..£il<^,„(£±l 

and the helicity-conserving amplitude 



ie] (120) 
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2x< 1 



t 



dx 



H{x,^,t) 
x-i 



dx 



2 



x-^ 
x + C 



2e 



2e 



Li, 



x + ^ 



+ l€\ -Li2(l) 



(121) 



The last expression presents the main result of this study. 

The transverse momentum squared, \P±\'^ > 0, can be 
expressed in terms of kinematic invariants as (j20p 



where 



(122) 



(123) 



is the minimal kinematically allowed value of the momen- 
tum transfer t. Note that the target mass corrections 
cx rri^ only enter our results through the dependence on 
^min and are always overcompensated by the finite-^ cor- 
rections (X |<| > |imin|- By this reason target mass cor- 
rections alone do not have any physical significance. 

The definition of the skewedness parameter ^ accepted 
in this work, Eq. (|18p . may not be the most natural one 
from the phenomenological point of view. Let 



Xb 



2-XB ' 



(124) 



where xb = Q^/{2pq) is the Bjorken scaling variable. A 
simple calculation yields for the relation of "our" f and 



t 



C= 1 + 7^ 



Q^n + iBt/Q-' 



(125) 



We assume that the GPD H{x,^,t) is continuous at 
the points x = ±^ while its first derivative can be sin- 
gular. The renormalization-group analysis suggests that 
the singularity in the derivative is at least logarithmic 

dxH{x, ^, t) ^ ln(x =F C) (if one does not employ some 
fine tuning) . It is easy to see that under this assumption 
the a;- integrals in Eqs. (|120p . (ll2ip converge and define a 
smooth function of ^ on the interval (0,1). This state- 
ment is obvious for the logarithmic terms, whereas in 
order to see this for the dilogarithm terms (which have a 
branch point at x = ^), it is convenient to rewrite these 
contributions as follows: 



1 



x~^ 



Li2 ( 



le] -Li2(l) 



x-S. Jo 



2« 



1-t 



■Ini ^ ln(i^ -zej.(126) 



For completeness we also give here an explicit expres- 
sion for the imaginary part of the helicity-conserving am- 
plitude: 



H{x,tt) 



2 



{^^+2\P^\^ed^)edi J dxH{x,^,t) 



1 1 



1 



2 x + ^ X - ^ 



In 



x + i 
2e 



(127) 



where we used that H{—x,£,,t) = —H{x,^,t) in order to 
bring the answer to this form. Note that the imaginary 
part depends on the GPD in the DGLAP region only, 
X > ^, in agreement with dispersion relations. 



V. NUMERICAL ESTIMATES AND 
DISCUSSION 

A detailed study of the numerical impact of kinematic 
power corrections goes beyond the tasks of this paper. 
This study has to be done at the level of cross sections, 
taking into account finite-t and target mass effects to 
kinematic (e.g. phase space) factors [20| and including 
the interference with the Bethe-Heitler process. Since 
the corrections to helicity amplitudes are calculated in 
this work for the simplified model case of a scalar target 
only, such a complete analysis is probably not warranted 
at this stage. 

In this Section we present numerical estimates for the 
kinematic power corrections to the imaginary part of the 
helicity-conserving DVCS amplitude, Eq. (|127l) . which 
gives the largest contribution to the cross section. To this 
end we use a model for the GPD H(x, t) corresponding 
to the A^ = 1 ansatz from Ref. [3l|. It is based on the 
so-called single-DD description which is defined by the 
"gauge" fixing condition 

(^f(J3,a,t) = I3g{l3,a,t) , 

imposed on the DDs / and g in Eq. (|59p . see Ref. [Slj 
for more details. It is assumed that the DD / takes a 
factorized form 



(128) 
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Here q{x, t 
we take as 



and 



0) is a (quark) parton distribution which 



3„Bt 



HP, a) 



3(1 
4 



(1 



(129) 



(130) 



The function /i(/3, a) satisfies the normalization condition 
/j'^i+^^l da h{f}, a) = 1. Note that we use q{x) ^ (1 — xY 
which is a standard approximation for the proton target, 
because this is case that is interesting phenomenologi- 
cally. For the pion one expects q{x) ~ (1 — x)^~^. 

We assume that the DD (|128p depends on t through 
the corresponding dependence of the Regge trajectory 

a{t) = 0.48 + 0.9 GeV^t 

and in addition involves a multiphcative factor e^* (see 
e.g. Ref. [13). 

The imaginary part of the amphtudc (jl27p involves 
H{x, f , t) in the region a; > ^ only. In this region one 
obtains a compact expression [3l| 



Hix,^,t)^>^ 



3x 



dp 



2n 



„Bt 



where Pi = [x — C)/(l — C) and P2 - 
We consider the following ratios: 

Im^(o) - lmA^°\ t 



(131) 
(x + 0/(1 + 0- 



ItoA 



(0) 
LO 



^-^ct{^,t) + ^c^2{^,t), (132) 



where 



(133) 



is the leading-order leading- twist result. Note that the 
e^* factor in the GPD cancels between the numerator and 
the denominator so that it is irrelevant for our purposes. 
The ratios Ct^m^{S,,t) still depend on t, however, because 
of the non-factorizable i-dependence of the GPD through 
the Regge trajectory. 

In Fig. [T] we show ct (0 1) a-s a function of ^ in the inter- 
val 0.001 < £, < 0.3 for f — (which also corresponds to 
a factorizable t dependence of the GPD) and for the val- 
ues —t = 0.2,0.5, 1.0 and 1.5 GeV^. The same is shown 
in Fig. [2] for the coefficient c,„2(0i). 

One sees that in the whole kinematic range |ct(0 01 ^ 
|cm2 (C; ^)l so that the target mass correction to the imag- 
inary part of the amplitude is almost negligible. The 
finite-t correction is, on the other hand, quite sizable. 
One has to have in mind that magnitude of the finite-i 
correction depends on the definition of the skewedness 
parameter. Using (|124p as a kinematic variable in the 
leading-order expression results in a redefinition of the 
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FIG. 1: The coefficient Ct{^,t) in the range 0.001 < C < 
0.3 for different values of t. The curves correspond to —t = 
(0.0, 0.2, 0.5, 1.0, 1.5) GeV^ 
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FIG. 2: The coefficient c^2{^,t) in the range 0.001 < C < 
0.3 for different values of t. The curves correspond to ~t = 
(0.0, 0.2, 0.5, 1.0, 1.5) GeV^. 



kinematic correction, which acquires a stronger ^ depen- 
dence but does not become smaller. 

For the particular case of the minimal allowed value of 
the momentum transfer t — tmin we find 



Imyl(") - ImA 



(0) 
LO 



ImA 



(0) 
LO 



-(0.62-0.65) 



t ■ 



(134) 

which holds with a good accuracy in a broad interval of 
^ and m^. 

To summarize, in this paper we have presented the first 
complete calculation of kinematic power corrections to 
the helicity amplitudes of deeply-virtual Compton scat- 
tering to the twist-four accuracy for a study case of 
a (pseudo)scalar target. Our main result is that both 



finite-t, ~ t/Q , and target mass. 



TO' 



7Q^, twist-four 



kinematic power corrections turn out to be factorizable, 
at least to the leading order in the strong coupling. Our 
numerical estimates using a certain simple model of the 
generalized parton distribution suggests that target mass 
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contributions to helicity amplitudes are very small, but 
the finite-i effects, on the contrary, rather significant in 
all kinematic regions of interest. The similar calculation 
for the DVCS from a nucleon is in progress and the re- 
sults will be published separately. 
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Appendix A: Flavor— singlet operators 

The generalization of the technique developed in 
Refs. [li [1] to the case of flavor-singlet operators is 
straightforward, so that we will not go into technical de- 
tails and discuss the key points only. Our notation in 
this Appendix follows Ref. [2|. 

The flavor-singlet quark operator of the leading twist 
0'^{zi,Z2) = 0+j^{zi,Z2) defined in Eq. pi)) has positive 
parity under charge conjugation and gets mixed by the 
renormalization group equations with the twist-two gluon 
operator 



{z,,Z2) ^fl+{zi7i)r+4z2n) . 



(A.l) 



The other existing quark- antiquark twist— 2 operator 
Ov{zi,Z2) + Ov{z2,zi)-Oa{zi,Z2) + Oa{z2,zi) has neg- 
ative C-parity. It does not mix with gluon operators and 
it also does not contribute to the OPE. The divergence 
of the flavor-singlet quark and gluon conformal opera- 
tors can be calculated following the procedure described 
in Ref. [2] . The result can be presented in the form of a 
scalar product 



{dOf 
(dO) 



N 

N ^ 



9 ' 



(A.2) 



where ~^q,~^g are the flavor-singlet (anti)quark-gluon 
and three-gluon nonquasipartonic operators, respectively, 
and 5* 1 2 J ^12 ^.re the corresponding coefficient func- 
tions. As discussed in detail in Ref. [1], the contributions 
of quasipartonic, alias four-particle twist-four operators 
can be omitted. 



The quark-gluon operators appearing in Eq. (jA.2[) 
present a flavor-singlet analog of the nonquasipartonic 
operators with "good" conformal properties introduced 



Q2{Z) 



v;+(zi)rv-(z3)+x-(^3rx+(^i)]/^+(^2), 
i?_+v;+(zi)r^_(z3) 

+ x-{zz)t^D^+x+{zi)\ r++{z2) ■ (A.3) 



Explicit expressions for the three gluon operators, ~^ , 
will not be relevant for the further discussion. 

The scalar product in (jA.2[) is determined by the re- 
quirement of hermiticity of the evolution equations. Its 
existence is ensured by conformal invariance. We need to 
know an explicit form of the scalar product in the quark- 
gluon sector only, and it turns out to be exactly the same 
as in the nonsinglet case. 

It is clear that {dO)\^ does not have a three-gluon com- 



ponent, so that '^\= Q. The wave functions ^ have the 

same form as in the nonsinglet case, whereas require 
a new calculation. The results read: 









JQ 




2 


W12 








Jo 



da aa {w2 — ^"3) 



a \N-1 



I a {w2 



da aa {w2 — w^^) 



a \N-2 



(A.4) 

-4 {N + 2) 



where = z/;2, W3) etc., 

and = 2 (iV - 1){N + 2){N + 4). 

Using orthogonality of the wave functions of different 
existing multiplicatively renormalizable operators with 
respect to our scalar product, cf. Q, the contribution 
of the divergence of the flavor-singlet twist-two operators 
{dO)%^ to the expansion of nonlocal operators Q , Q can 
be found as 



N 



^Jv|(||*2||'*l-*2(*2|*l))(50) 



where = (^\, ~^k)^ ^ = 1' ^ and the coefficient ^at is 



given by 



A 



N 



l*l|| 11*2 



|(*l|*2 



(A.6) 



Note that in order to calculate = and 

(v|/^|\p2) = it is sufficient to know the scalar 

product in the quark-gluon sector. 

The three quark-gluon operators in Eq. (jA.3[) are not 
independent (see Ref. 0) and their contribution to the 
OPE of the time-ordered product of two electromagnetic 
currents can be expressed in terms of Q2 alone. This 
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contribution, in turn, can be written in terms of a certain 
integral over the gluon position on the hght cone: 



dw {w - Z2)Q2{zi,W, Z2) ■ 



(A.7) 



Exphcit calculation shows that 



dw{w - Z2){^\)2{zi,W,Z2) = Cm\\^\ 



Z2 
Zl 



dw{w - Z2)i'^l)2{zi,W,Z2) = CAr(*i|^'2) - 



where 7^ and 7^ are the corresponding entries in the 
matrix of anomalous dimensions of the flavor-singlet 
twist-two operators. It follows then from (jA.5[) that the 
gluon operator (90);^ does not contribute to the OPE of 
two electromagnetic currents to our accuracy, while the 
flavor singlet quark operator (90)^ enters the OPE with 
the same coefficient as in the nonsinglet case. 



Appendix B: Translation invariance of the 
7?^— operator 

The formal proof of translation invariance of our re- 
sult [ll, m for the T-product of two electromagnetic cur- 
rents to the twist-four accuracy, Eq. pT|) . on the operator 
level is straightforward but rather long and technical. In 
this paper we have chosen to write this result in terms 
of the Ti{zi,Z2) operator defined in Eq. (1551) which is 
translation- invariant itself, Eq. (j57p . In this Appendix 
we give a proof of this relation, which holds up to twist- 
six corrections. 

Using the definition in Eq. ([55)1 expression on the l.h.s. 
of (1571) can be written in the form 



dw 



Z\2 



dw' w' — Z2 



Z12 Jz2 ^12 Zl 



-Z{w,w'), 



(B.l) 



where 



Z{w,w') = S-y-S+0i{w,w')-{Sa-l)02{w,w') 
= i^S+S- + So)Oi{w,w') - SoS.02{w,w') (B.2) 

and — dw + dw' ■ The operators ©1,2 are defined in 
Eq. ([S5]) and 0*-_^^ in Eq. (|^. Taking into account the 
following identity for the leading- twist projector 11 (see 
Ref. i): 

1 - 

n(2;,A)A"A" =x""n(a;,A) - -x^a"" / dull{ux,\), 

2 Jo 

(B.3) 

one finds 

S_0'^^{w,w') = n(a;,A)[i(Pn),0++(w,w')] 
= [i{Vx),0'=+{wy)] 
1 

- -x^ I duu02{uw,uw') . (B.4) 



With the help of this relation one easily obtains 

S-Oi{w,w') = [i{Px),Oiiw,w')]+ ... , 
S-02iw,w') = [i{Px),02{w,w')]+0i{w,w') 

duu [z(Pa;), 02{uw, uw')] 

+ ..., (B.5) 

where the ellipses stand for the contributions of twist- 
six operators. Finally, inserting these expressions in 
Z{w,w') and taking into account that 

So [ duu[i{Px),02{uw,uw')] = [i{Px),02{w,w')] 
Jo 

(B.6) 

one easily verifies Eq. ([57l) . 
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